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Abstract. A perturbed black hole emits gravitational radiation, usually termed the ringdown 
signal, whose frequency and time-constant depends on the mass and spin of the black hole. 
I investigate the case of a binary black hole merger resulting from two initially non-spinning 
black holes of various mass ratios, in quasi-circular orbits. The observed ringdown signal will be 
determined, among other things, by the black hole's spin-axis orientation with respect to Earth, 
its sky position and polarization angle - parameters which can take any values in a particular 
observation. I have carried out a statistical analysis of the effect of these variables, focusing on 
detection and measurement of the multimode ringdown signals using the reformulated European 
LISA mission. Next Gravitational- Wave Observatory, NGO, the third generation ground-based 
observatory, Einstein Telescope and the advanced era detector, aLIGO. To the extent possible 
I have discussed the effect of these results on plausible event rates, as well as astrophysical 
implications concerning the formation and growth of supermassive and intermediate mass black 
holes. 



1. Introduction 

Astrophysical observations to date have provided sturdy evidence that black holes (BHs) may 
exist and play an important role in many physical processes [H [2l |3]. With direct evidence 
still lacking, it is expected that observation of gravitational waves (GWs) from merged BHs will 
not only provide indisputable evidence for the existence of BHs, but also the ability to extract 
accurate information about the progenitor system and the BH. 

The ringdown radiation consists of a superposition of, in principle, an infinite number of 
essentially damped sinusoids, termed quasi-normal modes. Their frequencies and time-constants 
depend only on the mass and spin of the BH - a consequence of the no-hair theorem [3]. In 
a recent work [5| we have argued that the amplitude terms of the various quasi-normal modes 
encode important information about the origin of the perturbation that caused them, such as 
the component masses of the progenitor binary. This allows performing parameter estimation 
on the system from the strong-field regime, as opposed to using the inspiral phase. However, in 
that study, as well as in previous studies of parameter estimation from ringdown signals [4, 6j, 
only a small region of the parameter space was explored. In a realistic scenario, a BH ringdown 
signal can have any sky location and polarization. All the while, it could have originated from 
a BH with any spin-axis orientation with respect to Earth. These variables have a significant 



effect on the observed ringdown signal and a direct impact on the science we can achieve by 
observing GWs from merged BHs. 

In the present study, I have investigated in detail how these angular parameters affect the 
detection and measurement of the ringdown signals. To this end, I have varied the angular 
parameters over their full range, thereby considering a large population of BH mergers. I have 
considered supermassive BHs (SMBHs of mass > IO^Mq) visible in NGO and intermediate mass 
BHs (IMBHs of mass ~ IO^Mq) observable in ET and aLIGO. I have computed the probability 
distribution functions of signal to noise ratios (SNR), as well as measurement errors of a chosen 
set of parameters, for a wide range of the BH mass and for mass ratios between 2 and 20. 
Finally, I translate these probabilities to proportions of observed events in NGO and ET that 
will yield parameter errors below certain thresholds and discuss how observations of ringdown 
signals could help in dealing with open questions on the existence and history of SMBHs and 
IMBHs. 




Figure 1: Noise power spectral densities for various 
proposed configurations of LISA-like space detectors. 
The dashed red line corresponds to the original LISA 
mission, while the rest of them refer to European 
designs of LISA. In all cases, the galactic binary 
white dwarf confusion noise [3 [S] is included, which 
has a negligible visible effect on the newer LISA 
curves though, due to their, almost two orders of 
magnitude worse sensitivity. Additionally, a low 
frequency cut-off - not shown - was induced at 
5 X 10~^iJz. In this study, I am using the latest 
arrangement for the European mission of LISA, 
labeled eLISA-NGO in this graph. 



Figure 2: This plot shows the relative amplitudes of 
modes (3, 3), (2, 1) and (4, 4) over (2, 2) as a function 
of the mass ratio. The diamond points show the 
values that correspond to a time lOM after the peak 
luminosity of the 22 mode in the equal mass case, 
(see also Fig. 3 in [S]) Fits to these points were used 
in this previous work. On the other hand, the circle 
points, which are used in the present study, were 
computed by taking into account all the points in 
the waveform, in a time region starting at the peak 
luminosity of 22 and ending 30M later. The solid 
lines shown, constitute fits to these circle points, 
given by expressions (|6|)-(|8|. 



2. Full population Analysis 

I have closely followed the procedure of Ref. [5] for estimating the signal-to-noise ratios and 
measurement errors. In this Section I will discuss the signal model and the parameter space 
covered in this study. 

In the generic case where we have a network of detectors, we write the response to a ringdown 
signal as: 

/i"(t)= i?Lne-*/"'-cos(u;,„„t + 7Ln)> (1) 

£,m,n>0 



where the superscript a is an index denoting the detector in question and uJimm T£mn ai'e 
the frequencies and time-constants of each mode, which are functions of the mass and spin 
magnitude of the BH. For further reference, see [9l [10]. In this study I neglect modes with 
overtone index n > 1, thereby considering only the least damped modes. From now on, the n 
index (n=0) is omitted. The terms and ^-re the following combinations of the antenna 
pattern functions -F^, F^, amplitude factors aim and the angular functions Yf^{i), Y^"^{l) : 

BL = ^\l{F'^yln" + {FiY^n\ (2) 
itm = 4'em + m(l) + tan 

Here, (pim are arbitrary constant phases of each mode. The effective amplitudes B^m vary 
inversely with the luminosity distance and proportionally to the intrinsic amplitudes a£m of the 
modes, which are determined by the numerical simulations. The angular functions (i) are 
the following combinations of the spin- weighted spherical harmonics [B]: 

Y^"'ii) = _2l^''"(^o)-(-l/_2r'-'"(^o). (4) 

The antenna pattern functions are functions of the sky location coordinates, 9 and (p a-nd 
the polarization angle tp, that is, F^{6, ip, ip), F^{9, ip, ip). The spheroidal harmonics are angular 
functions of the inclination angle, l and the azimuth angle (p. The first refers to the angle formed 
by the BH's spin angular momentum and the line-of-sight, while the latter is the azimuth angle 
defined in a non-rotating frame fixed to the BH. 

2.1. Chosen waveform 

The ringdown waveform used is of the form described by Eqs. ([T])-([3]). It is a signal comprised of 
four modes, with mode indices {i, m, n) = (2, 2, 0), (3, 3, 0), (2, 1, 0), (4, 4, 0). Our choice is based 
on the ordering of the various modes according to power output, as determined from numerical 
simulations of non-spinning unequal mass binaries [21 [TTl [3] . 

Numerical simulations of merging black-hole binaries were performed using the BAM code 
[m [13], so as to obtain the amplitudes, of the various modes in Eq.Q, as well as their 
frequencies and time-constants (see also Table I of [5]). These simulations involve the case of 
initially non-spinning BHs in quasi-circular orbits and for different mass ratios of the binary. 
For the mass ratios, q = {1,2,3,4} the simulation results were first presented in \14:\ [T5] I16j . 
while an additional simulation of a g = 11 binary was carried out in [5]. 

The amplitude terms in Eq. [2| are given by the expressions: 



a22{q) = 0.25e-'?/^•^ (5) 

a33(g) = 0.18022(9) (g-l)°-'', (6) 

a2iiq) = 0.15 022(9) (g-l)°■^^ (7) 

044(9) = 0.018022(9) 9°-^^ (8) 



These constitute fits, that were produced by fitting the merger- ringdown part of the numerical 
simulations data, taking into account all the different mass ratios for which these were performed. 
All points in a time region beginning at the peak luminosity of 22 and ending 30M later were 
considered, for each mass ratio. As opposed to the method that was applied in our previous 
work [5], where the relative amplitude values at lOM were used, this approach is expected to 
be more robust and to average out any numerical noise that might be present in this part of 
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Figure 3: Frequency distributions for the total SNR and the measurement errors in NGO, for a 5 x lO^M© 
BH situated at a luminosity distance of 6.73 Gpc {z ~ 1). The top left plot shows the probability to 
have a detector-BH configuration which will yield the SNR shown, while the rest of the plots concern 
cumulative frequency distributions for the measurement errors. In each graph, the comparison is shown 
among different mass ratios, q, specifically taking the values of 2, 10 and 20. For q=2 - black solid 
lines - the parameter estimation for the BH mass and spin is outstanding in all configurations. This 
effectiveness degrades considerably with increasing q. For q of around 10 it is still acceptable, while at 
q~20, all parameters except the mass are very likely to have huge errors, of the order of 200%. 



the waveform. Nevertheless, as can be seen from Fig. [2j these two methods do not give very 
different results. Note also that the above fitting functions, as well as the mode frequencies 
and time-constants, may be less accurate in the higher mass ratio values of around 20, where 
extrapolation has been performed. 

2.2. The simulations 

In Ref. [5] we had ignored the effect of the various angles {9, ip, ip, l} on the quasi-normal 
mode spectrum and their impact on the detection and measurement of ringdown signals. To 
assess this effect, the aforementioned analysis was repeated by varying the angular parameters 
{6, 99, i}- Specifically, six uniformly spaced values were chosen for these angles. This results 
in 6^ = 1296 distinct relative orientations between the detector and the BH and its spin axis. 
Additionally, a couple of simulations with eight uniformly spaced values were performed, that 
is 8^ = 4096 configurations, to allow the comparison. Hence, it was decided that six values in 
each parameter was acceptable in capturing the behaviour of the observable quantities. 

The values in the polarization angle, ij: and the azimuth sky location angle, ip were linearly 
sampled in the ranges [0, vr] and [0, 27r] respectively. Whereas, the values in the inclination angle 
i and sky position, 6 - which range from to vr - are deduced from the uniformly spaced values 
of cos(t) and cos{6) in the range [— 1,+1]. Note that this excludes configurations of optimally 
oriented binaries, that is of i = and /, = vr. 

2.3. The parameter space 

The parameter set of the ringdown signal in the case of a non-spinning black hole binary, 
consists of the following nine parameters: {M, j, (7, 9^, 4>}- Namely, the mass M, 

the dimensionless Kerr parameter or spin magnitude, j of the BH, the mass ratio q of the 
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Figure 4-' Same as in Fig. ^ but concerning a 25 x IO^Mq BH. This figure concerns a BH 5 times more 
massive than that of Fig. ISTbut some of the results are actually slightly worse, as the multimodal signal's 
power spectrum is shifted away from the lowest part of NGO's sensitivity curve. 



progenitor binary, the sky location vector (Dl, 0, ^) of the BH with respect to Earth, the 
polarization angle ip, the inclination angle l and the BH azimuth angle (p. Note however, that in 
this case, the final spin of the BH is directly mapped to the mass ratio of the progenitor binary, 
therefore q and j are not treated as independent. 

The above parameters are the standard ones, as they pertain to all of the modes. Additional 
parameters can be introduced, that are characteristic to each mode, such as an initial phase 
factor (pem, see Eq. ([s]). Therefore, the total number of parameters can increase with the number 
of modes. We are considering a four mode signal, therefore we have a total of 13 parameters. 
One thing to note is that the i = m modes have a nearly consistent rotational phasing, while 
the i ^ m modes seem to have somewhat distinct associated dynamics, with differentiated 
amplitude and phasing during the merger process [llj . 

By virtue of the large number of parameters involved, it was unmanageable to treat the effect 
of all of them in this analysis. Thus, some of the above mentioned parameters had to be fixed. 
Specifically, the initial phase angles, (p£mO in all the four modes considered, were plainly chosen 
as zero. In addition, the luminosity distance is chosen to be 6.73 Gpc for NGO, and 1 Gpc for 
ET and aLIGO. Lastly, the BH azimuth angle, (p was given the value 7r/3. Note that this angle 
does not have an effect on the SNR, but of course needs to be considered in the waveform. 

2.4- Choice of masses and mass ratios 

The BH mass and the binary mass ratio constitute key parameters and the results depend quite 
strongly on them. The reasoning behind the choice of mass values is the following. First of all, 
the low and high mass end is limited by the sensitivity band of the detectors. The existence of 
BHs in the mass range 10^ — IO^Mq, is highly predicted by the mass - velocity dispersion in 
the galactic bulge of low-mass, low-luminosity galaxies, as well as in a number of galaxies which 
contain active galactic nuclei \17\ I18|. The evidence for SMBHs ranging from IO^Mq to 10^ Mq, 
is quite abundant. They are thought to dwell at the centers of most galaxies. For various recent 
SMBHs mass estimation results and methods see for instance \19\ [20l [2H [221 E3 • The SNR in 
NGO is quite low for BH masses of less than about 5 x IO^Mq. Therefore, I take the lower end 
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Figure 5: Same as in Figs. |3]and|4j but for a IO^Mq BH. The results slightly deteriorate with respect to 
the lower mass BH of Fig. ffl 



of the mass range to be 5 x IO^Mq and consider two other values of 25 x IO^Mq and a lO^M© 
BH, to cover the interesting range of masses potentially observable in NGO. 

In the case of ET and aLIGO I have considered three IMBHs of mass 2OOM0, 6OOM0 and 
1OOOM0. It is believed that BHs in this range are situated in the centers of many globular 
clusters. However, their existence is being questioned, the evidence is thought to be strong but 
circumstantial [Ml O ISSl [261 [27] . 

Concerning the mass ratios, only unequal mass binaries are presented, except for one case 
in aLIGO. We do not consider the equal mass case as such systems are not as likely to occur 
in nature as slightly asymmetric ones. For NGO we examine the mass ratios of 2, 10 and 20, 
whilst lower mass ratios of 2, 5 and 10 are considered for ET. We, therefore, have in total, 9 
different sets of simulations for each detector. Let me emphasize here that there is a possibility 
that massive BH binaries in the early universe {z < 10) will have a mass ratio significantly 
larger than one. For instance, in [28] it is suggested that low-redshift massive BH mergers occur 
predominantly with a mass ratio of 10 or higher. 



3. Signal detectability and Parameter estimation 

I will discuss the total signal to noise ratio, as well as the fractional errors in estimating the 
following 5 parameters: {M, j, q, Di^, l}. These errors are actually the quantities, ax = \/C\x, 
which are computed from the covariance matrix, Cim \29[ [30l [311 [32j . 

The results from all the distinct arrangements of the system (see Section 2.2) are presented 
via cumulative frequency distribution plots. That is to say, the different system configurations 
are classified according to the value they render in the error of the observable quantity in 
question. The proportion of a number of occurrences in a small width of values should, to a 
good approximation, equal the probability that a randomly placed observer will measure that 
quantity to take this range of values. 
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Figure 6: Similar graphs as those in Figs, [s] - [sj although for the ET detector, concerning a IMBH of 
2OOM0. The luminosity distance of the BH is now picked closer, at 1 Gpc {z ~ 0.2). Also, note that the 
progenitor mass ratios now take lower values, namely 2, 5 and 10. 



3.1. NGO 

The sensitivity curv ej^that we use is what is thought to be possible for NGO and it is contrasted 
with other LISA-like sensitivity curves in Fig. [T} It refers to a 4-link interferometer, comprised of 
one mother and two daughter spacecraft, having armlengths of 10^ m and trailing a few degrees 
behind the Earth, in heliocentric orbit. 

The main noise contributions are the acceleration noise, the shot noise, as well as some other 
measurement noise. These are respectively: 

SaccAf) = 1.37 X 10-32(1 + — ^2^^-l, 

SsN,m = 5.25 X 10-23 rn^Hz-^, 
SoMN,m = 6.28 X 10- 23 H z'^ . 

The formula for the amplitude sensitivity curve is, 

4S^) = V5-|t(/)^^1^^^±^^±^ Hz~-I\ (9) 

while the transfer function is 



T{f) 



i+IttttWtI- (10) 



with L = lO^m and c = 299, 792, 458 metres per second. 
3.2. Results for NGO 

Our results for NGO are plotted in Figs, [s p] For the 5 x IO^Mq, 25 x IO^Mq and IO^Mq BHs 
the SNR curves (top left subplots of Figs. |3j[5j) peak at around 300, 600 and 1700 respectively. 

^ see https:/ /lisa-light. aei.mpg.de/bin/view/DetectorConfigurations 




Relative frequency distributions were chosen for their plotting, as they portray clearly where the 
maximum occurs, as well as how they fall off. The higher mass ratio SNR curves have a similar 
outline. The fact that all of the curves resemble log- normal distributions, with steep risings and 
long tails, is mostly attributed to the fact that only a small fraction of the configurations, those 
close to the optimal orientation of the binary, will yield the highest SNRs. 

Regarding the rest of the plots in Figs. [3]j5| the general trend is that the mass and spin have 
comparable, low fractional errors, whilst the other group of parameters, namely the luminosity 
distance, mass ratio and inclination angle, yield an order of magnitude higher error values. 
This is not surprising, considering that the mass and spin have a direct effect on the modes' 
frequencies, w^mn and time-constants, T£mm quantities that determine to first order the shape 
of the ringdown waveform. 

Quoting a few numbers for a progenitor of mass ratio 2, the probability to get a binary merger 
event that will yield a mass error less than 1% is correspondingly 58%, 55% and 10%, as we go 
from the lower mass to the higher mass value. The lowest mass value actually gives the best 
results, whereas the signal's main power content takes place near the lowest sensitivity area of 
the NGO detector. The spin magnitude exhibits a similar trend with the probabilities to fall 
below a 10% fractional error being 96%, 90% and 67% respectively. 

For the second group of parameters, I again quote how likely it is to do better than 10%. 
For the luminosity distance, the values are 38%, 41% and 20%, corresponding to the BH masses 
5 X IO^Mq, 25 X IO^Mq and 10^ Mq. The mass ratio is harder to determine accurately, with 
the likelihoods being 32%, 15% and only 2% respectively. Lastly, there is a 40%, 37% and 13% 
likelihood of achieving an accuracy better than 10% in the inclination angle. 

3.3. Einstein Telescope and Advanced LIGO 

ET's very low sensitivity curve accounts for impressive results in the mass range ~ 5OOM0 to 
~ IOOOMq and for mass ratios between 1 and 5. I consider the sensitivity curve designated 
ET-B [33J, whose noise power spectral density is given by Sh{f) = 10~^''/in(/)^ Hz~^, with: 

hnif) = 2.39 X 10-27 349 ^-2.145 

+ 1.76x-°-^2^ 0.409x1-^°, (11) 



where x = //lOOHz. 

As for advanced LIGO, the noise spectral density is 



10l6-4(/-7.9)2^ 08^-4.69 

Hz-\ (12) 

where x = //215Hz. 

34. Results for ET and aLIGO 

The results obtained for ET and advanced LIGO are plotted in Figs. [6]|9j We fix the luminosity 
distance of the BH to be 1 Gpc. For the lowest mass considered (a 200Mq BH) although the 
SNR could be pretty high (in the range 30-200), errors in the estimation of parameters are poor 
(see Fig. [6]). {DL,q,L} have 50% probability to be measured to an accuracy of ~ 50%, while 
errors on {M,j} are 90% and 60% likely to be below 10%. For the higher mass ratios the results, 
as expected, are worse. 

The results for 6OOM0 and IOOOMq BHs are shown in Figs. [7] and |8j The parameter 
estimation accuracies for these systems observed with ET is almost as good as that for a SMBH 
with NGO. Referring to the heaviest BH and mass ratio in the range 2-5, it is 100%-99% and 
96%-53% likely to acquire errors below 10% for the BH mass and spin. For the {DL,q,L}, the 
efficiencies are correspondingly 50%-15%, 23%-0% and 60%-17%. The results are similar for the 
6OOM0 case (see Fig. [7]). 

Fig. [9] shows two examples for aLIGO: a BH resulting from an equal mass binary and one for 
which the mass ratio is 5. Most configurations give ringdown SNR values in the range 10 to 30, 
not large enough for a good estimation of parameters. This translates to a 30% likelihood for 
a fair measurement of 10% accuracy in the mass and mass ratio, while the luminosity distance, 
spin and inclination angles are all measured to accuracies far worse than 10%. 

4. Astrophysical implications 
4-1- Supermassive black holes 

The presence of SMBHs in the centers of massive galaxies seems to be a well established fact. 
Detecting their gravitational wave signals can give additional clues on their spatiaj^ and mass 
distribution, as well as help discriminate among the different scenarios of their formation and 
growth. This could be done, for instance, by measuring the BH mass function as a function of 
red-shift. Additionally, determining the mass ratios [5] of these early universe merger events will 
be an important piece of information in selecting out the current models on SMBH formation. 
Note that, in NGO, BHs of mass higher than 10^~*Mq are visible almost entirely due to the 
ringdown signal that they emit rather than the inspiral signal. 

Several studies have been realized in the field of predicting the coalescence rates of SMBHs 
[3ll[35l[36ll371[38l[39lllQlllI]. Let us admit an event rate of ~ 10 yr^^ at z ~ 1, which ascertains 
a scenario of most efficient BH coalescence [42j. Then, assuming that the BH masses are both 
around 5 x IO^Mq, Figs. [3]and[4| there is a good chance that in 6 of the events, the BH mass 
will be measured more accurately than 1% and that in 9 of the events the BH spin magnitude 
will have an accuracy better than 10%. As for the parameters {D^, q, l}, in approximately 3-4 
of the events they will feature errors lower than 10%, while in 7 of the events, the errors will lie 
below 30% for and below 20% for q and i. 

^ This fit was provided by C. Capano, Syracuse University and is tuned for detecting binary neutron stars. 

^ Unless the inspiral phase is inside NGO's band, it will not be possible to determine the sky position of the 

signal. The only hope in this case would be the existence of an electromagnetic counterpart to the merger event. 



Shif) = 10- 



+ 123.35 



1 - 0.23x^ + 0.0764x4 
1 + 0.17x2 



0.2 
0.18 
0.16 
0.14 
0.12 

0.1 
0.08 
0.06 
0.04 
0.02 
















— 


q = 2 


















q = 5 


















q = 10 










































































\ 





































200 400 600 800 1000 1200 1400 
total SNR 



1 

0.9 
0.8 
0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 














/ / 






/ 










/ / 






J- 


1 - q = 2 
q = 5 



0.1 1 
fractional error in D, 



1 

0.9 
0.8 
0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 


































f / - q = 2 




/ ■■■■ q = 5 




/' q = 10[ 


0.01 0.1 



fractional error in M 




1 

0.9 
0.8 
0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 
























f f 












J ^ / 
















- q = 2 




^ J 




q = 5 




jj 




q = 10 



1 

0.9 
0.8 
0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 










































J / 








- q = 2 


~i 1 1 1 1 iTi 




■■■■ q = 5 




q = 10 



0.01 0.1 

fractional error in j 



0.1 1 

fractional error in q 



0.1 1 
fractional error in i 



Figure 8: As in Figs. [6]and[7| but for a lOOOM© BH, a sweet spot in ET. The results are very encouraging 
for this size of BH and distance, with a very high probabiUty to get errors fower than 10% in most of the 
parameters at mass ratios 2-5. 



Let me clarify that ringdown signals in NGO, will be detectable out to about z ~ 5, but the 
distance of z = 1 was chosen to make the quote on the errors. If we take the distance at z = 3, 
an optimistic rate of merger events will be of the order of 100 [32]. If we assume that they all 
involve BHs of masses ~ 2-3 x lO^M©, then in around 90 of them the erroij^in the mass will be 
below 10% and in 60 of them the spin magnitude error will fall below the 10% threshold. For 
the rest of the parameters, about 50 of the observed events should yield measurements better 
than 50% accuracy. Although this looks poor at the outset, it should suffice for a statistical test 
of diff'erent models of BH formation and growth. 
4-2. Intermediate-mass black holes 

The existence of intermediate-mass black holes remains uncertain, as is their mass distribution. 
Colliding globular clusters in interacting galaxies could be a mechanism to obtain a compact 
binary IMBH system |43l E]. Another possibility could be the formation of a binary IMBH 
inside a young dense stellar cluster, especially when the fraction of binary stars is adequately 
high ||26|. 

Estimates of IMBH-IMBH coalescence rates can be found in |441 [26l O US]. A relatively 
optimistic rate is Ropt = 0.007 GC~^Gyr~^ , where it has been assumed that 10% of star clusters 
are sufficiently massive and have a sufficient stellar binary fraction to form an IMBH-IMBH 
binary once in their lifetime, taken at 13.8 Gyr. The maximum possible rate would come from 
assuming that all of the star clusters satisfied the above conditions. The corresponding rate 
value is then Rmax = 0.07 GC~'^Gyr-^. 

If the number of relatively luminous galaxies within a distance of 1 Gpc is approximately 
5.3 X 10^, [36] and the number of young dense stellar clusters per such galaxy is of the order of 
100, then the optimistic estimate gives 0.037 events per year, while the maximum rate would be 
0.37 events per year. An event within this distance will in some likelihood, involve IMBHs with 
masses between ~ 6 x lO^M© and ~ W^Mq which means a relatively fair chance for errors to 
be low for several of the parameters using Einstein Telescope (see Figs, ulpl). 



An additional simulation at z = 3, of a 5 x 10® BH, with a progenitor mass ratio of 2 was performed. 
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Figure 9: Frequency distributions involving the most optimistic scenario for advanced LIGO: a IOOOMq 
BH which is the merger of an equal mass binary, that is g = 1. The BH is again situated at 1 Gpc. A 
mass ratio of 5 is shown as well for contrasting. 



5. Conclusions 

The present study constitutes a sensible and realistic approach to the subject of parameter 
estimation from a multimodal ringdown signal, inasmuch as it is supposed to be emitted from 
a merged binary in a generic configuration. Parameters such as the inclination angle l, along 
with the sky location 6 and if and the signal polarization ■0 have an effect on the observed 
quasi-normal mode spectrum. Their impact on the detectability and parameter estimation has 
been assessed, by performing a large sample of Fisher-matrix analysis simulations, allowing for 
a simple statistical analysis of the results. 

I am quoting frequency distributions for the errors at the representative distances of z ~ 1 
(6.73 Gpc) for supermassive BHs, in NGO, and of z ~ 0.2 (1 Gpc) for intermediate-mass BHs, 
in the Einstein Telescope and advanced LIGO. An additional simulation was performed for a 
5 X 10® Mq BH at z ~ 3. The results are quite satisfactory in determining the mass and spin, 
especially for the low mass ratios from 1 to ~ 5, where in typically 90% of the cases their 
errors fall below 10%. The effectiveness in measuring the luminosity distance, mass ratio and 
inclination angle is almost an order of magnitude worse. 

The effects of these results on, as much as possible, realistic event rates in NGO and ET are 
discussed. The likelihood to have a waveform parameter measured to an accuracy of a certain 
threshold translates to the same proportion of observed events featuring error values below that 
threshold. As an example, if supermassive BHs within a luminosity distance of z = 1 coalesce 
at a rate of ~ 10, then NGO could act as a SMBH dynamics probe, as almost all of these events 
will yield very low errors 1% — 10% in the BH mass and spin, while in half of the events the 
errors in the luminosity distance, mass ratio and inclination angle will be of the order of 20%. 
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